Nucleon correlations and Higher Twist effects in nuclear structure functions 
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The overlap of the nucleons in nuclei plays an important role in understanding the nuclear depen- 
dence of deep inelastic scattering data. It is shown that the nuclear modification of the higher twist 
scale can be essentially determined by the overlapping volume per nucleon and this effect gives a 
large contribution to nuclear shadowing for small x and low Q^. In this kinematical region there is 
also a moderate enhancement of the longitudinal structure function in nuclei. 



o 
o 

o 

Q 



> 
in 

o 

O 



X 



The interpretation of the HERA and RHIC data and 
the future experiments at LHC require a deep under- 
standing of the behaviour of the parton distributions in 
the small x region both for the free nucleon and for nuclei. 

The main ideas of the dynamical description of the 
small X physics are the saturation of the parton den- 
sities and the shadowing effects. Both these problems 
can be studied from the fundamental point of view by 
the non linear QCD evolution equations |H. However, 
since they are related with confinement, some comple- 
mentary approaches, as classical gluon field theory j2| or 
the pomerons exchange , take into account the non 
perturbative dynamics. 

Nuclear shadowing is a crucial ingredient of the initial 
state effects in hard phenomena in proton-nucleus and 
nucleus-nucleus collisions and then some useful , model 
independent, phenomenological analyses of the nuclear 
modifications of the parton distributions have been pro- 
posed 4J. 

The major difficulty for a clear theoretical description 
of the nuclear shadowing is that the available data on 
deep inelastic scattering (DIS) on nuclear target in the 
small X region are also for small and then the standard 
QCD evolution becomes unreliable. 

Recently the calculation of the nuclear parton distri- 
bution functions in the Gribov approach has been up- 
dated and the failure of the leading twist model to 
describe the DIS data indicates significant higher twist 
(HT) contributions in the small x region. In particu- 
lar, the HT contribution to nuclear shadowing turns out 
about 50% of the whole suppression and the authors ar- 
gue that in the kinematical range covered by data it is 
possible to have this large corrections due to diffractive 
effects. On the other hand, in ref. the HT contribu- 
tion turns out small ( see later). 

In this letter I suggest that in the small x and low 
region the HT contributions are relevant for the agree- 
ment with the present day data on nuclear shadowing. 
This conclusion, that can be also obtained by resum- 
ming QCD power corrections 0, is based on an intuitive, 
coarse-grain analysis of the HT effects by the overlap of 
nucleons in the nucleus. 

When nucleons overlap, the effective average volume 



per nucleon decreases due to nuclear density and correla- 
tions. This reduction of the average interquark distances 
, due to overlapping volume per nucleon V^^ in the nu- 
cleus A, is an effective way to describe a greater number 
of quarks in a larger bag with , however, a reduced ef- 
fective volume per nucleon, Va, with respect to the free 
nucleon volume V^v- 

This idea has been apphed to the EMC eflFect ,by the 
Qf-rescaling 9] model, and also to the nuclear shadowing 

I shall show that by a simple analysis, based on the 
overlapping volume per nucleon in the nucleus, one can 
suggest that: 

i) there is a nuclear modification of the HT scale which 
gives a relevant contribution to nuclear shadowing in the 
present day data; 

ii) the nuclear correction to the longitudinal structure 
function, f j"^ , are not so large as preliminary reported 
in ref. ^. 

Let us first discuss the higher twist effects which are 
subasymptotic scaling violations, 0{y? /Q'^), whose typ- 
ical scale fj? is fixed by the matrix elements of the rel- 
evant operators in the light-cone expansion of the prod- 
uct of the currents operators. In particular, Jaffe and Sol- 
date 01 found a complete set of nine operators, totally 
symmetric, traceless in Lorentz indices and depending 
on color, spin and flavor indices, which permits us to 
evaluate the twist four, spin two, contribution to electro- 
production in QCD. 

The HT formalism has been largely applied ^3 and 
it requires a specific model to evaluate the non pertur- 
bative, target dependent, matrix elements. For example, 
in the MIT bag model the twist four correction to the 
nucleon structure function F2N turns out 
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where and Vn are respectively the nucleon mass 
and volume, ac is an effective coupling constant in the 
bag ( usually ac — 0.5 is assumed) and Ki — 1.08 and 
K2 = 0.17 are integrals over Vm of the bagmodel spinors 
for massless quarks in the lowest modes |l2j . 

This result, by itself, is not useful for the comparison 
with DIS data because one needs to analyze the x de- 
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pendence and, moreover, in the small x, low region 
the valence approximation of the bag model is not re- 
liable. A complete calculations of the matrix elements is 
still missing and, on the other hand, it is quite useful to 
have, at least, a phenomenological understanding of the 
HT contributions, crucial for the agreement with data at 
small Q^. 

Then, following the structure of eq.(l), let us write 
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and try to have some indication of the effective weight of 
the matrix elements for determining the HT scale of the 
nucleon in the small x region. The next step will be to 
generalize the analysis for a nuclear target. 

However, before entering in the phenomenological de- 
tails, let us consider some results which follow from gen- 
eral arguments. 

Let us write the F2 structure functions as a leading 
twist term (LT) plus a Oil/Q'^) HT contribution. For the 
nucleon case we shall approximate F2N = F2D /2 since at 
high energy the nuclear effects in deuterium are negligible 
and from now on the formulas refere to structure func- 
tions per nucleon. For deuterium one has F2D — F2^ + 
F^j^ and for a larger nucleus A, F2A = F^J^ + F^/. By 
defining R^^ = F^J/F^^ and Ra = F2A/F2D ,to order 
0{1/Q^), one obtains 
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Then the HT corrections can suppress the (LT) term only 
if RaF^j^ > F2A -Since in the small x region Ra < 1, 
if the HT correction increases the (LT) term (i.e. it is 
positive) then the effective HT scale must be larger in 
deuterium with respect to the larger nucleus A . Vicev- 
ersa, if the HT term are negative the larger nuclei have 
a larger effective scale at small x. 

Concerning the sign of the HT correction to F2D let us 
notice that: 1) In the MIT bag model this overall correc- 
tion is negative ( see eq. (1)); 2) in the parametrization 
of the deuteron structure function data in ref. the 
0{1/Q^) leading correction in the small x region is neg- 
ative. 

By assuming this indication, one concludes that the 
nuclear HT corrections are negative and there is a larger 
HT effective scale at small x with respect to the nu- 
cleon case. Moreover one can use the phenomelogical in- 
put given in ref. [l^ and write that 
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where ci = l.bmGev^ and ca = -g.SSSGew^ [iJ.Of 
course this approximation is valid only in the small x 
region considered in the data of ref. |lj| . 



Let us now consider the nuclear effects. The HT cor- 
rections come essentially from the matrix elements of 
quark-quark-gluon and/or four quark operators [l^ in 
the target and then in a nucleus there are peculiar con- 
tributions, with respect to the free nucleon, due to the 
nuclear binding and to partons coming from different nu- 
cleons in the nucleus. A rigorous calculation is quite 
hard, but, according to the parametrization in eq.(2), 
it is natural a description of the nuclear effects which 
takes into account the nuclear binding energy by an ef- 
fective mass of the nucleon, M* = Mm + e (where e is the 
removal energy [Tsf.and the correlation among nucleons 
by their overlap in the nucleus as previously discussed. 
Within this extremely simplified approach of the nuclear 
effects, the HT contribution to F2A per nucleon can be 
parametrized as ( see eq.(2)) 
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where Va/Vn = 1 — /Vn and Ca{x), related to the 
matrix element of the relevant operator in the nucleus, 
takes also into account the HT contribution due to par- 
tons from different nucleons.lt is reasonable to assume 
that this dynamical effect is proportional to the overlap- 
ping volume per nucleon, , and then to write 



Ca{x)=Cd{x)+C{x) 



Vn 
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where C{x) is independent on A. 

Then the HT nuclear correction turns out that 

F2A -F2D JJ^^[l+lix)—]=F2D^A 



(7) 



where 7(0;) = Ca{x)/Cn{x) will be approximated with a 
constant parameter, 7, in such a way that the HT nuclear 
scale is completely determined by the overlapping vol- 
ume per nucleon. By combining the previous eqs.(3,7), 
the relation between R^ and Ra in terms of the HT 
corrections is given by 
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which is the starting point for the comparison with the 
experimental data. 

In the previous formula the only free parameter is 7. 
Indeed the overlapping volume per nucleon has been eval- 
uated in ref. by including a Reid-soft core potential 
between nucleons and considering only the two nucleon 
overlap. In ref. 'I6] the previous evaluation has been 
reanalized with a reduction of a factor ~ 0.65 with re- 
spect to the initial values. Then the overlapping volume 
per nucleon for different nuclei is fixed as the 65% [l^l 
of the values obtained in ref. 0. The removal energy e 
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FIG. 1: Fig.l Comparison of the leading twist theory plus 
higher twist contributions according to eq.(8) to the NMC 
data [l^on Fic jF-iu^ sees text. The statistical and systematic 
errors are quadratically added. 



is small ( few Mev) and negative and it has a negligible 
quantitative effect. F^j^ is given in eq.(4) and F2B is ex- 
perimentally known for small value of x and accurately 
parametrized 14]. Finally, one has to know BJ^ and dif- 
ferent LT inputs change the fitted value of the parameter 
7- 

From this point of view, it should be clear that the 
agreement with data depends on the relative weight be- 
tween the LT and the HT contributions. The power cor- 
rections are an important element of the whole analy- 
sis but their magnitude is correlated with the the initial 
parametrization of the structure functions for the evo- 
lution. In ref. 1^ a model of diffraction is used in order 
to obtain an initial condition for DGLAP evolution at 
Qq — 4 Gev^. On the contrary, in ref. 7] a model, to 
all twist orders, is proposed for the full low region 
and there is no QCD evolution. In this latter framework, 
where a small HT correction is reported, is not simple to 
handle the 0{1/Q^) and clearly separate a "LT" term. 
Then, in the present analysis, one considers the lead- 
ing twist results in the Gribov approach obtained in ref. 

In Fig.l is reported the comparison of the leading 
twist theory plus higher twist contributions , according 
to eq.(8) to the NMC data ISjon i^2c/-F2D-The dashed 
and dotted-dashed lines are respectively for 7 = 8 and 
7 = 4 and the dotted line is the leading twist result in 
ref.Q. Fig. 2 is the same comparison with the Ca-40 data 
and in Fig. 3 the comparison is with the ratio F2pb/F2c- 
In fig. 4 has been considered the dependence at fixed 
X = 0.0125 of the ratio i^2Sn/-F2C |I3- 

Let us now comment on the longitudinal structure 
function In ref. 0| has been reported a quite 

large value ( about a factor 4) of the ratio between the 
longitudinal and transverse structure functions in N-14 
for X ~ 0.0045 and low with respect to the deuterium 
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FIG. 2: Fig. 2 Comparison of the leading twist theory plus 
higher twist contributions , according to eq.(8) to the NMC 
data [Hon F2C a/F2D, see text. The statistical and systematic 
errors are quadratically added. 
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FIG. 3: Fig. 2 Comparison of the leading twist theory plus 
higher twist contributions to the NMC data on F2Pb / F2C ,see 
text. The statistical and systematic errors are quadratically 
added. 




FIG. 4: Fig. 2 Comparison of the leading twist theory plus 
higher twist contributions to the NMC data on F2Sn/F2c at 
X = 0.0125 and different Q^,see text. The statistical and 
systematic errors are quadratically added. The dotted line is 
the leading twist result 0.975 0. 
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case. 

This ehancement has been investigated by many au- 
thors d m 113 with the conclusion that a large HT 
correction for nytrogen is required. Indeed, also in this 
case, the effect is due to a balance between LT term and 
HT correction which depends on the absolute value of 
the HT scale rather than on the ratio between the ni- 
trogen and deuterium HT scale. By the present analysis, 
an enhancement should be found because the HT nuclear 
scale increases. However it turns out that the effective 
HT nytrogen scale at x c± 0.0045 is much smaller than 
the the value of O.ObGGev'^ required to obtain the 
reported large correction. 

Moreover one has to remember that:l) the matrix ele- 
ments which produce the HT correction for the Fl and 
F2 structure functions are correlated but are not the same 
and then the previous results on the HT scale in nuclei 
cannot be directly applied to the longitudinal structure 
function; 2) the data in ref. |0| are under reanalysis [2^ . 

Let us conclude with few comments. 

The question if the observed nuclear shadowing is a 
leading twist or a higher twist effect is still open. 

In the present analysis, the HT contribution to nuclear 
shadowing is relevant, according to ref. 0, but this could 
depend on the initial parametrization. A conclusive an- 
swer requires an updating of the phenomenological fits 
which includes the 0{1/Q^) corrections, the QCD evolu- 
tion and also the experimental data on the longitudinal 
structure function. 

The overlap of the nuclcons is an important element 
in understanding the nuclear shadowing ( also at leading 
twist 10]) and the EMC effect 9]. This can be probably 
considered as a first signal towards a geometrical criti- 
cal transition in heavy ions collision at high energy and 
density, associated with percolation 23] rather than with 
the temperature. 
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CERN Theory Division for hospitality. 
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